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Abstract 

It is proved that the attenuation function of a viscoelastic medium 
with a non-decreasing and concave creep compliance is sublinear in the 
high frequency range. 
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Notation. 

K the set of real numbers 

C the set of complex numbers 

fog composition {f ° 9){x) f{g{x)) 

^[/] (P) = hp) Laplace transform of f{t) f {p) = /~ e-P* J{t) dt 

f{k) Fourier transform of /(a;) f(k) = J_^e^^'''^ f{x) dx 

1 Introduction. 

Ultrasonic tests of mechanical properties of polymers and bio-tissues assume 
that the material is viscoelastic. In applications a linear viscoelastic material 
is always assumed to have a completely monotonic (CM) relaxation modulus 
[U [5J 13] . This assumption has important consequences for the attenuation and 
dispersion of acoustic waves. In particular, the increase of attenuation with fre- 
quency in the high frequency range is slower than linear [4] . This fact cannot be 
verified experimentally because experimental test cover the range of frequencies 
below 250 MHz. The results of ultrasonic tests of polymers and bio-tissues sug- 
gest a power-law dependence of attenuation in the frequency range accessible 
to measurements, A{uj) = Auj°' with a > 1. This power-law dependence is 
then extrapolated to high frequency range, where it contradicts various general 
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principles. It is then explained in terms of ad hoc lossy wave equations [5J O [7] 
and wave equations with fractional Laplacians |151 116] which are incompatible 
with viscoelastic constitutive equations. 

Superlinear growth of attenuation in the high frequency range is inconsistent 
with finite speed of propagation. The requirement of finite speed of propaga- 
tion imposes an upper growth limit uj / \niuj) on high-frequency attenuation: the 
propagation speed is finite if w/ln"(a;), a > 1 but not for a < 1. Superlinear 
growth of attenuation in the high frequency range is also inconsistent with vis- 
coelastic constitutive equations even if the propagation speed is infinite, for 
example it is inconsistent with Newtonian viscosity, where the attenuation sat- 
isfies the power law const x w^/^, and with strongly singular locally integrable 
completely monotonic relaxation moduli [4 . 

According to the theory of wave propagation in viscoelastic media with com- 
pletely monotonic relaxation moduli, developed in [HIH], in the low frequency 
range the wave attenuation in a viscoelastic solid increases like W"'^^", < a < 1, 
while a viscoelastic fiuid is characterized by a rate of increase oj", < a < 1. 
Low frequency behavior of attenuation is accessible to measurements and is con- 
sistent with observations. In the high frequency range the growth of attenuation 
is however sublinear. High frequency attenuation is also sublinear in viscoelastic 
media with spatial derivative operators of fractional order |14j . 

Although practical modeling of viscoelastic media routinely assumes that the 
relaxation modulus is completely monotonic, there is hardly any a priori evi- 
dence for this assumption. The only general physical principle - the fluctuation- 
dissipation theorem - implies that the relaxation modulus is completely positive 
[9] . The consequences of complete positivity for the dependence of attenuation 
on frequency are however at present unknown. In this paper we explore the 
consequences of a different extension of the complete monotonicity assumption, 
based on the assumption that the creep compliance is non-decreasing and con- 
cave. Non-negative functions with these two properties are said to belong to 
the CrF class, which is an abbreviation of Priiss' unfortunate term "creep func- 
tions". The above assumption, proposed and explored in Jan Priiss' book |10] . 
is much weaker than the hypothesis of a completely monotonic relaxation mod- 
ulus, assumed in [SlU]. The latter is equivalent to the creep compliance being 
a Bernstein function, which is defined in terms of an infinite sequence of in- 
equalities for the derivatives of the creep compliance of arbitrary order. In the 
CrF class only the three first of these inequalities are required to hold. The 
Bernstein class and the CrF class includes Newtonian viscosity characterized by 
linear creep. Unlike the Bernstein property, the CrF property of a function can 
be easily visually recognized from its plot. Furthermore, the CrF property is 
stable with respect to perturbations which are small in the sup norm. 

It will be shown here that wave attenuation in media with CrF creep com- 
pliances is bounded by a function a -I- 6 A sharper form of this bound was 
obtained for creep compliances in the Bernstein class in [51 |3]. The theory 
developed in [8l H] relies on a spectral analysis of attenuation and dispersion. 
Unfortunately a spectral analysis of attenuation and dispersion similar to the 
completely monotonic case has not been possible but it is very likely that the 
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bound obtained here for the CrF class of creep comphances can be significantly 
sharpened. 

After a recapitulation of necessary facts from the theory of CM functions 
and Bernstein functions (Sec. [5]) the essentials of Priiss' theory are recalled in 
Sec. ini The main result concerning the bound on wave attenuation is derived 
in Sec. [51 

2 Some properties of completely monotonic and 
Bernstein functions. 

Definition 2.1. A real function f on ]0,oo[ is said to be completely monotonic 
(CM) if f £ C" and (-1)" D" f{x) > for for every integer n>0. 

Theorem 2.2. (Bernstein) 

A real function f on ]0, oo [ is CM if and only if there is a positive Radon measure 
H on [0, cx)[ such that 



A CM function is locally integrable if and only it is integrable over the 
interval [0,1]. The set of locally integrable CM (in short LICM) functions is 
denoted by £. 

Theorem 2.3. HJl \^ 

A real function f is LICM if and only if there is a positive Radon measure fj, 
such that 



and eq. ([T]) holds. 

Definition 2.4. A real function f defined on [0, oo[ is said to be a Bernstein 
function (BF) if it has a derivative which is a LICM function. 
The set of Bernstein functions is denoted by the symbol *8. 

It is obvious that a Bernstein function is non-negative, non-decreasing and 
concave. 

From Bernstein's theorem the following representation of Bernstein functions 
can be derived: 

Theorem 2.5. fT^ 

A function f : M+ — R is a Bernstein function if and only if the are two 
non-negative real numbers a, b and a positive Radon measure A on R+ such that 




t > 



(1) 




(2) 




(3) 
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and 

f{t)=a + bt+ [ [l-e-*^*] A(dr), t>0 (4) 
Inequality (jSj is equivalent to the following pair of inequalities 

rj/(dr) < oo, / i^(dr) < oo (5) 

04] "']l,cx)[ 

Theorem 2.6. /// satisfies eq. ^ with X satisfying eq. thenlimt^aa fit)/t = 
b. 

Proof. Expansion in power series shows that (1 + t)"^ /2 < e* < (1 + t)"^ for all 
t > 0. Hence (1 + i)"^ < e"* < 2(1 + i)"^ for t > and, upon integration, 
t/il + i) < 1 - C-* < 2</(l + t). 
Now 

— = - + 5 + / A(dr) (6) 

t ^ "']0,oo[ ^ 

The integrand of the integral on the right-hand side of eq. ^ tends to as 
< — )■ oo. The integrand is also bounded from above by 2r/(l + rt) and by 
2r/(l-|-r) for i > 1. In view of eq. ([3]) and the Lebesgue Dominated Convergence 
Theorem the integral on the right-hand side of eq. (|6]) tends to as i — cx). The 
thesis follows from eq. ([5]). □ 

Theorem 2.7. fW^ 

If the real function f defined on ]0,oo[ is CM and g € S, then the composition 
fog is CM. 

If f, g are Bernstein functions on [0, oo[, then the composition fog is a Bernstein 
function. 

If a positive real function f defined on [0, oo[ is a BF then 1// is CM on ]0, oo[. 



3 The CrF class of functions. 

In every neighborhood of a bounded CM function in the space of bounded 
continuous functions C°([0, cx)[) endowed with the sup norm there is a function g 
which is not CM, for example g{x) = f{x)+sm{x)/n. Hence the CM property of 
the relaxation modulus cannot be experimentally confirmed by a direct method 
although it can be inconsistent with experimental data in a specific case. Thus 
complete monotonicity has to be considered as an a priori constraint on the 
interpolation of experimental data by matching them against a multi-parameter 
family of completely monotonic functions (Prony sums, Cole-Cole or Havriliak- 
Negami functions, the Kohlrausch- Watts- Williams function). 

It is therefore interesting to note that the linear upper bound on the high 
frequency behavior of the attenuation function is not a consequence of complete 
monotonicity of relaxation moduli. As an alternative to viscoelasticity based 
on completely monotonic relaxation functions, a theory based on an experimen- 
tally verifiable concavity property of creep compliances will be presented here 
following [10 . 
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(a) 



Figure 1: (a): The function defined by ([7]), for a = 1/2 (dashed), a — 3/2 
(solid) and a = 5/2 (dot-dashed). (6) The 3rd order derivative D'^/ for the 
same values of a. 

Definition 3.1. Let I be a segment of the real line. 
A function / : / — > R is said to be concave if 

^f{x) + (1 - ^?) f{y) <f{^x + {l- ^) y) 

for all x,y e I, <■& <1. 

Bernstein functions are non-negative, non-decreasing and concave. 

Definition 3.2. A function f : M+ ^ M is said to be a CrF if it is non-negative, 
non- decreasing and concave. 

O. denotes the set of all the CrFs. 

Every Bernstein function is obviously a CrF. 

Example 3.3. The function 

fix) Pe-y^dy (7) 
Jo 

(Fig. op zs a CrF for a > 0, but it is not a BF if a > 1. Indeed, if a > then 
/, /' > while /" < 0. 

It a<l then f'{x) = e-="° is CM, hence / G *B. 

If a > 1 then f"'{x) = aa;"^^ [a a; — {a — 1)] e~^ changes sign at x = 
{a — l)/a. Hence f is not a Bernstein function in this case. 

Theorem 3.4. Every CrF f can be expressed in the form 

f{t) =a + bt+ [ k{s)ds (8) 
Jo 

where a, b are non-negative numbers and the locally integrable function k is non- 
negative, non-increasing and fc(s) — > for s — > oo. 
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Proof. Let a := /(O) and b := inft>o f{t)/t. The function b assumes finite values 
because concavity 



implies that 



f{dt) 1-d fit) 
-a > 



dt - t 

for t > and < i9 < 1, hence f{t)/t is decreasing and has a finite limit b. The 
function / is non-negative and non-decreasing, hence it is bounded on an interval 
[0,e], e > 0. Since — / is convex, / has a derivative /' at all but a countable 
set of points. The derivative /' is almost everywhere non-negative (because / is 
non-decreasing) and almost everywhere non- increasing (because / is concave). 
The second-order derivative /" exists and is non-positive almost everywhere. 
Hence f'{t) < f'{t) + t f"{t) = d[t f'{t)] and thus, upon integration and dividing 
by t, f'{t) > f{t)/t — a/t. Since f{t) is non-decreasing, the right-hand side of 
this inequality is > b. 

The function k[t) :— f'(t) — b is measurable, almost everywhere non- negative, 
and non-increasing. Hence k is locally integrable near (otherwise / would be 
infinite at 0). In view of monotonicity \mit^ock{t) exists and, in view of the 
definition of &, it vanishes. □ 

□ 

If / is a CrF function satisfying eq. ([8]) then D/ exists almost everywhere 
and D/ = k, /(O) = a and limt^oo D/(i) = b. 

Examples 3.5. 

1. The function f{t) := a+bt+ct" , a,b, c > 0, < a < I, is a CrF function, 
with k(t) ^at"-'^. 

2. The function f defined above is not a CrF if a > 1 . 
Theorem 3.6. 

(i) If f e Q and g{x) f{x) then g e *8. 

(ii) If f CI ^ then there is a CrF g such that f{x) = x'^ g{x). 
Proof Ad (i) 

Eq. ([8]) implies that / has the Laplace transform 

f{p) = - + — + -Hp) 

p p'^ p 

hence 

Hp) = ap + b + pk{p) 

But 

P kip) =p e-P'' fc(r) dr^ d [l - e'P''] fc(r) = - / [l - e'^'P] dfc(r) 
Jo Jo Jo 
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Define the Radon measure A on R+ by the formula / A(dr) = — / Lp{r) dk{r) 
for continuous functions ip with compact support in R_|_. The measure A is 
positive and 

-A(drj ~ — — = / — —dr < oo 



]o,oo[ ^ + r y]o,oo[ 1 + r Jo (1 + ^) 

We have used the facts that hmr^o+ rk{r) =0 (because k is locally integrable) 
and k is non-increasing. Consequently the Radon measure A satisfies inequality 
([3]). We now have the identity 

/•oo 

pk{p) = / [1 -e-''P] dfc(r) 

We thus see that the function g has the integral representation ^ and is thus 
a Bernstein function, q.e.d. 

If the limit k^ := limr-i.o+ k(r) exists then u can be defined as the Borel 
measure i^(]r, s]) = k{r) — k{s) for < r < s, provided k is redefined as a 
right-continuous function. 

Ad (ii) 

The function / satisfies eq. ^ with the Radon measure A satisfying inequality 
(O. Since r/{l + r) > 1/2 for r > 1, it follows from eq. ^ that for every r > 
the right-continuous function k{r) := A(]r, oo[) — is non-negative 

and finite. Since 

b pb 

k{r)dr= / r X{dr) + b k{b) - a k{a) 



is in view of eq. ([5]) i finite for < a < 6, the function k is locally integrable on 
It is also non- negative, non- increasing and, on account of ([3]), it tends to 
zero at infinity. The integral in eq. ^ can be integrated by parts: 



(1 - e-P'') dfc(r) = p / e-f k{r) dr = p k{p) 

0,oo[ "']0,oo[ 

Hence f{t) = p^ g{p) with g{t) ■= at + b + k{s) ds and g G Q by Theorem[ 

□ 

Remark 3.7. The definition of a complete Bernstein function (Appcndix\S^ is 
equivalent to the statement that, assuming again the equation g{x) = f{x), 
g ^ d "if and only if f <E'^. 



4 The wave number function in materials with 
a CrF creep comphance. 

We now consider a viscoelastic material whose creep compliance J is a CrF. 



7 



Consider the Fourier and Laplace transform of the viscoelastic equation of 
motion assume 

Pp'^{l= ~pG{p)Pu (9) 

and note that pG{p) = l/[pJ{p)] [llj . Define the wave number function 

p G C, in such a way that k = —iK{p) and Re k{p) > is a solution of eq. dH). 

In view of eq. ^ the wave number function can be expressed in terms of 
the Laplace transform J of the creep compliance 

K{p)^p'/^p'/^ [p^J{p)]'^' (10) 

If J S Q then, by Theorem 13. 6[ f{p) := p^ J{p) is a Bernstein function. Hence 
K ^ [pf{p)Y/'^ e /Cq := 251/2^ where 03" := {/" | / S 25}. We note that 
*B" C *B for < a < 1 because /" is a composition of two Bernstein functions 
a: — >■ and /, see Theorem 12.71 However pf{p), f G ®, is in general not a 
Bernstein function, for example if /(p) = 1 — e~^. 

In |3] it was proved that /C<8 S" H p^^^^ is the set of all wave number 
functions compatible with the assumption that the relaxation modulus is LICM, 
or, equivalently, that the creep compliance is a Bernstein function [llj . The 
set of Bernstein functions *B is a subset of £3. We thus expect that the set 
JC<3 C /Cq. In order to check this fact directly we shall need the following 
theorem (Proposition 7.11 in [13 ): 

Theorem 4.1. If < a < 1 then ^ {g e ^ \ p^'"" g{p) e d}- 
Corollary 4.2. 5- n p^/^ 5^ = p^/^ 

Proof. Theorem |X1 and Theorem |X3] imply that p^/^ ^ n p^/^ gr 

If / e STlp^/^ ^ then f{p) = p^/'^ g{p) with g e^. Theorem lO implies that 
g e hence / e p^/^ ^ g ^j^ □ 

□ 

It follows from Corollary 1121 that d f] p^^'^ ^ = p^^^S^^^ C p^/^ S^/^ as 
expected. 

Returning to the case of a creep compliance in the CrF class, note that the 
wavenumber function has the form k{p) = p^^^ /(p)^^^, where / has the integral 
representation 

f{p) = a + bp+ [ [1 - e-f '■j y[Ar) (11) 

J\0,oo[ 

with a, 5 > and a non- negative Radon measure v satisfying the inequalities 

Using Theorem l2.6l it is now easy to prove that k{p) — h^l^ p + i?(p), where 
i?(p) = o[p] for p CO. The phase function in Green's functions has the form 

exp(]9< - K(p)|a;|) = exp(p(i - 5^/^ |2:|) - R{p)\x\) 

with p — — icj, w G M. The arguments used in the proof of Theorem 7.2 in [4] 
can be used to demonstrate that the Green's function vanishes for |a;| > 
hence is the wavefront speed. 
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5 Main theorem. 



Behavior of the attenuation and dispersion function for real frequencies (imagi- 
nary p) is of particular interest. We shall therefore examine the behavior of the 
analytic continuation of a Bernstein function f{p) on the imaginary axis. 

Analytic continuation of / to the imaginary axis yields a function F{uj) := 
/(— icj), w e R. /R,(aj) := Re F{—iuj) and /i(w) := Im F{—iuj). Hence 

fn{uj) = a + / [1 - cos(ra;)] i^^dr) + [I- cos(ra;)] ^(dr) (12) 



/l(^) 



sin(rti;) iy(dr) + 



sin(ra;) v^dr) 



(13) 



Lemma 5.1. The integrals in eqs. I11S^13\} are bounded by linear functions of 
the variable tu. 

Proof. Note that | sin(a;)/x| < 1 (the numerator and the denominator vanish at 
and the absolute value of the derivative of the numerator does not exceed the 
value of the derivative of the numerator equal to 1). Hence | sin(rti;)/r| < |a;|. 
In the same way one shows that |(1 — cos(rw)/r| < Hence 



sin(ra;) i^(dr) 



0,1] 



sin(raj) 



0,1] 



rj/(dr) 



< \oj\ I r ^{dr) 

10,1] 



and 



]o,i] 
Furthermore 



[1 — cos(ra;)] ^(dr) 



1 — cos(ra;) 



0,1] 



riy{dr) 



< \u)\ I rj/(dr) 
]o,i] 



sin(rw) j/(dr) 



[1 — cos{tuj)] i^(dr) 



< I iy{dr) < oo 

]l,oo[ 

< 2 / i^(dr) < oo 

']l,oo[ 



Hence |/R(a;)| < a + C\uj\ and |/i(w)| < 6 + | < D\uj\ for some positive reals 
a,b,C,D. □ 

Theorem 5.2. 



\ki{uj)\ <K + L\lu\ 
|kr(w)| <K + L\ij\ 



(14) 
(15) 



for some real constants K and L. 
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Proof. If TO = maxja, C}, M ~ niax{6, D} then 

|F(c.)| = V/r(w)' + < 72 (to + Af 

and therefore 

\K{-ioj)\ < 2i/4 (to + M < 2^/4 (a^ + Ml/2 1^1^ 

where TV := to/(2M1/2). The thesis foUows with K = 2^/^ TV and i = 21/* Afi/^. 

□ 

It foUows that the attenuation is majorized by a hnear function of frequency. 
In the high frequency range it increases at most at a hnear rate. This bound 
can perhaps be sharpened. 

6 Concluding remarks. 

It has been shown that the assumption that the creep comphance is non- 
decreasing and concave imphes that the attenuation function increases at a 
sublinear rate in the high-frequency range. This result is somewhat weaker 
than for materials with completely monotonic relaxation moduli and Bernstein 
class creep compliances [8j |4] . 

This result shows that some acoustic wave equations applied to ultrasonic 
tests of materials and especially for ultrasound in bio-tissues still popular in 
literature [5j |6l [7] are not compatible with viscoelasticity. Such equations are 
based on the assumption of a power law attenuation with an exponent exceed- 
ing 1. We have thus shown that such equations are incompatible with the 
assumption that the material is linear viscoelastic and the creep compliance is 
non-decreasing and concave. 

The results obtained in this paper are compatible with experimental data. 
In the frequency range covered by measurements the frequency dependence of 
attenuation is represented by a power law with an exponent 1 < a < 2 over less 
than three decades of frequency expressed in MHz. This fact is consistent with 
the results obtained in [4] for low-frequency attenuation in viscoelastic solids. 
It should however be borne in mind that the sound absorption measured in 
heterogeneous media includes loss due to scattering which is here not accounted 
for. 
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A 



Definition A.l. A real function f on [0, oo[ is said to be a complete Bernstein 
function if there is a Bernstein function g such that f{x) — g{x). 

f{x) = a + bx + x [ (16) 
^]o,oo[ x + r 

Theorem A. 2. ^ 

Every CBF can be continued to an analytic function F on C_ satisfying the 
following conditions: 



(i) F assumes real values on . 



(a) a finite limit limx-s-o F{x) exists and is non-negative; 
(Hi) sgn(Ini [F{z)]) — sgn(Im [z]); 

(iv) The function F has an unique integral representation 

f tz — 1 , ^ 
F{z) ^a + bz+ / p{dt) 



0,oc[ 



t+ Z 



with a,b > and a positive Radon measure p satisfying the inequality 
satisfying the inequality 



j 

J]0,oo[ 



< QO 

l+t 



If a complex valued function F satisfies either the conditions (i)-(iii) or 
condition (iv) then the restriction of F to K+ is a CBF. 

This theorem has the fohowing corohary: 

Theorem A.3. /// is a CBF and < a < 1 then f{x)°' is also a CBF. 

Theorem A. 4. // /, g G and < a < 1 then the pointwise product h := 

Proof. The functions f,g map M+ into R+, hence h has the same property. 

The functions f,g are holomorphic in C± and map C± into C±, hence h{z) 
is holomorphic in C+ and < arg h{z) < a tt + (1 — a) tt = tt. By Theorem IA.2I 

heS. □ 
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